Euler Coefficients and Restricted Dyck Paths 

Heinrich Niederhausen and Shaun Sullivan, 
Department of Mathematics, Florida Atlantic University, Boca Raton 



We consider the problem of enumerating Dyck paths staying weakly 
above the x-axis with a limit to the number of consecutive \ steps, or 
a limit to the number of consecutive /" steps. We use Finite Operator 
Calculus to obtain formulas for the number of all such paths reaching a 
given point in the first quadrant. All our results are based on the Eulerian 
coefficients. 

1 Introduction 



One of the most recent papers on patterns occurring k times in Dyck paths 
was written by A. Sapounakis, I. Tasoulas, P. Tsikouras, Counting strings in 
Dyck paths, 2007, to appear in Discrete Mathematics [5]. The authors find 
generating functions for all 16 patterns generated by combinations of four 
up (y) and down (\) steps. A Dyck path starts at (0, 0), takes only up and 
down steps, and ends at (2n, 0), staying weakly above the cc-axis. Returning 
to the X-axis at the end of the path has the advantage that every path 
containing the pattern uduu, say, k times, will contain the reversed pattern 
ddud also k times when read backwards. This reduces significantly the 
number of patterns under consideration. Dyck paths containing k strings 
of length 3 were discussed by E. Deutsch in [1]. 

In this paper we consider only the patterns and d^ , for all integers r > 
2, and we will investigate only the case k = 0, which means pattern avoid- 
ance. It has been shown in [5] that the generating function f{t) for avoiding 
x''(or d^) satisfies the equation / {t) = l+Y^lZ^ t^f [tf — — 



-t-rf{tr 



l-2t 



How- 
ever, we will allow the Dyck paths to end at (n, m), m > 0, which removes 
the above mentioned symmetry, as shown in the following two tables. 
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The number of Dyck paths avoiding uuuu 
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The number of Dyck paths avoiding dddd 

The two tables indicate the differences between the two problems, both 
starting out from equal counts on the x-axis (m = 0). Because only points 
(n, m) with n + m = mod 2 can be reached by a Dyck path, we consider 
the lattice points (2n + 6, 2m + fe), for b — 0,1. We first show that the 
number of Dyck paths to (2n + b, 2m + b) avoiding d^ equals 



Dyck {2n + b, 2m + b; d'' 



2m + 5+1 / n + m + b 
n + m + b+ l\ n — m 



1 



where the Euler coefficient [5] is denoted by 



n + m + + 1 
n ~ m 



^ ^_j^^Jn + m + 5 



1=0 



m + 5 - 




(2n- 


Vb- 


ri\ 


i 




\n — 


m — 


ri 1 



(see Definition [4] and expansion (fl^ ). More about Euler coefficients can be 
found in Section HI For given m, the number of Dyck paths 
Dyck (n + m,m ^ n; d^) to (n + m, to — n) avoiding d^ has the generating 
function (over n) (1 - T)" (1 - t)"™"^ {rV (1 - i) - (1 - T) {2t - 1)), as 
shown in (fTT|) . Note that the coefficient of t™ in this generating function 
equals the original Dyck (2m, 0; d^). 

Next we show that the number of Dyck paths to {2n + b, 2m + b) avoid- 
ing u'^ equals 



D (2m + b,2m + b;u'^) 
1 



2m+&-l 

y 

^ n + TO + 6+ l 

1=0 



- 2m - b 
I \ i 



n + m + b+ l— i 
n + m + b — i 



except for the original Dyck path counts to (2n, 0), which either must be 
gotten from those to {2n — 1,1), or from the Dyck paths to (2n, 0) avoiding 
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(F . The case r — A seems to be very special. We conjecture in Section IXTl 
that in this case the generating function for the Dyck paths equals 



Dyck (4m — n — 1, 2m — n + 1; u**) 
= (^3 + t - V(l + i)'+4i3) (^) " /2, 



hence Dyck 




Throughout the following sections we will discuss ballot paths (weakly 
above y = x), with steps | and —>■, instead of Dyck paths. The trans- 
formations D (n, to) = Dyck (n + m,m — n) and Dyck {2n + h, 2m + 5) = 
D (n — m, n + m + b) , with D {n, m) counting ballot path to (n, to), go back 
and forth between the two equivalent setups. Of course, the pattern be- 
comes the pattern f"^, or iV, and (T becomes or . 



2 Ballot paths without the pattern 

Deflnition 1 s„(m;r) — Snim) is the number of {1,^} paths staying 
weakly above the diagonal y = x from (0, 0) to (n, to) G avoiding a 
sequence of r > consecutive —^ steps. We get sq (to) ~ 1 for all to > 0. 
We set Snim) — if n < orifm+l~n>0. 

Lemma 2 T/ie following recurrence holds for all m > n > 0: 

Snim) = s„_i(to) + s„(m - 1) - Sn-r{m - 1). (1) 

Proof: The number of paths reaching (n, m) is obtained by adding the 
number of paths reaching [n — 1, to) and (n, to — 1), but subtracting paths 
that would have exactly r — > steps. Those forbidden steps occur necessarily 
at the end of the path, so they are preceded by an up step, and must come 
from (n — r, TO — 1). ■ 

We now extend s„(to) to all integers m by first setting so(m) — 1 and 
using ID) to define the remaining Sn{m) for m < n — 1. 

Lemma 3 (s„) is a polynomial sequence with degs„ = n. 

Proof: We proceed by induction on n. Clearly, deg (sq) = 0. Suppose 
Sk{m) is a polynomial of degree k for < < Then s;+i(m) — si+i{m — 
1) = si(rn) — s;_r+i(m — 1), which implies the first difference of si^i{m) is 
a polynomial of degree I. Thus, S(_|_i(to) is a polynomial of degree / -f 1. ■ 
By interpolation we can define (s„) on all real numbers. 
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Definition 4 The Eulerian Coefficient is defined as 



(see lil^) ). Note that for r — 2 the Euler coefficient equals the binomial 



coefficient 




The following table shows the polynomial extension of s„ (m) . The num- 
ber of {t, — >} paths to (n, m) avoiding a sequence of 4 ^ steps appear above 
the y = X diagonal. The numbers on the diagonal (n, n), 1, 1, 2, 5, 13, 36, . . . , 
are the number of Dyek paths to (2n,0). Of course, s„ (n) < C„, the rt-th 
Catalan number. 
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The path counts s„ (m) and their polynomial extension (r — 4) 



Theorem 5 

X — n + lfx + 1 



Snix) = 



x + 1 



[n/rj 



X — n + 1 v-^ , ^ f X + 1\ f X + n — ri 



Proof: We saw that (s„ [x)) is a basis for the vector space of polynomials. 
Using operators on polynomials, we can write the recurrence relation as 



= B-B''E-^ 



(2) 



where B and are defined by linear extension of Bsn{x) = Sn-i{x) and 
E°-Snix) = Sn{x + a), the shift by a. The operators V = 1 - E~^ and E~^ 
both have power series expansions in D, the derivative operator. Hence 
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B must have such an expansion too, and therefore commutes with V and 
E"-. The power series for B must be of order 1, because B reduces degrees 
by 1. Such hnear operators are called delta operators. The basic sequence 
(6„ {^))n>o ^ delta operator _B is a sequence of polynomials such that 
deg hn = n, Bbn (x) = (x) (like the Sheffer sequence s„ (x) for B), and 
initial conditions 6„ (0) = (5o,nfor all n £ Nq. In our special case, the basic 
sequence is easily determined. Solving for in ([2|) shows that 

E^^J2^'- 

i=Q 

Finite Operator Calculus tells us that if 1 + (j{B), where a{t) is a 

power series of order 1 [31 (2.5)], then the basic sequence bn{x) of B has 
the generating function 

^bn{x)t" ^{l + a{t)r. 
n>0 

Thus, in our case bn{x) = [i"] {1+t + t'^ ^ h f-^Y = (^^ . Since the 

Sheffer sequence (s„) has initial values s„(n — 1) = Sn^Q, using Abelization 
[3] gives us 

, , X — n + , , x-n + l/x + l\ 
X + 1 X + 1 \ n / ^ 



Corollary 6 The number of Dyck paths to (2n, 0) avoiding r down steps 

IS 

3 Ballot paths without the pattern 1'^ 

Definition 7 tn(rn;r) = tn{m) is the number of {T,^} paths staying 
weakly above the line y = x from (0, 0) to (n, m) avoiding a sequence of 
r > consecutive ] steps. We set t„(m) = Q if n < Q orm+l=n>Q. 

This time we do not immediately have a polynomial sequence, as the 
table below shows. The path iV-i [EN'-'^)'' to (fc, (r - 1) (fc + 1)) is the 
only admissible path reaching the point (fc, (r — 1) (fc + 1)) (all others would 
have r or more iV-steps). Hence tn-i {{r ~l)n) = 1 for all n > 1, and 
tn-i (m) — for m > {r — 1) n. The only other I's in the table occur in 
column 0, to (jn) — 1 for to = 0, . . . , r — 1, and for all other values of m. 



5 



The table contains a strip weakly above the diagonal y — x where 

tn {m) = tn{m - I) + tn~i {m) (5) 

This happens for 0<n<TO<7i + r because paths in this strip cannot 
have r consecutive vertical steps. All paths that reach a point (n, ni) for 
m > n + r and violate the condition of not containing must have this 
pattern exactly at the end of the path, which means that they end in the 
pattern EN"^ . Hence for m > n + r we get the recurrence 

tn {m) = tn{m - 1) + t„-i (m) - t^-i {m - r) (6) 

We assume that t„ (m) = for all to < n (also for n 

We can find a recursion that holds for all to > n as follows: For n>\ 
we always have tn (n) = (n), because i„ {n — 1) — 0. From ([5]) follows 
by induction (inside the exceptional strip) that tn (m) = J2iLn ^n-i («) for 
all rt < TO < n + r. For the values of m on the boundary of the strip we 
have 

n+r n+r 

t„ (n + r) = ^ i„_i (i) - i„_i (n) = ^ (?) 

i—n i—n+1 

from ([5]) and ([6]), and after that by induction using ([6]), 

tn (w) — X]i=m+i-r ^"-1 (*) ^'^^ — n+r. We cau write both recursions 

together as 

m 

tn {m) = ^ i„_i (i) (7) 

z— max{n,m+l — r} 

for all m > n. We can avoid the difficulty with the lower bound in the 
summation by setting t„ (m) = for all m < n. Call the modified numbers 
t'^ (to). The new table follows the recursion 

m 

for all TO > rt > 0. The 'lost' value t„ (n) can be easily recovered, because 
tn (n) tn-i (n) = 4-1 ('^)- 
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The restricted ballot path counts t„ (m) (r = 4). 



In order to show the polynomial structure in the above table, we transform 
it into the table below by a 90° counterclockwise turn, and shifting the 
top I's flush against the y-axis. In formulas, we define Pn{m) = ^ 
l)m — n) for to (r — 1) > n > (or (to) = p(r-i)(n+i)-m ("- + !))■ The 
recursion t'^ (to) X]"=m+i-r ^'i-i (*) 'along the previous column' becomes 
now a recursion p„(m) = X][=o Pn-i (™ ~ 1) 'a-long the previous row'. More 
precisely, for (r — 1)to — 7i>77i— 1 > 1 , i.e. , to > 2 and n < (r — 2) to 
holds 

(r — l)m — n 

p„(m) = 4^_i((r- l)m-n) = ^ C-2 (*) 

■i— (r — 1) (m— 1) — n 
(r — l)m— n r — 1 

= X! P{r~l){m^l)-t{m-l) =^Pn-]{ni-l){8) 

i=(i — l)(m — 1) — n i=0 

The numbers Pn (to) for < n < (r — 2) m are exactly the cases where 
(m) is positive, and the only additional numbers needed in the recursion 
([5]) are the numbers 

P(r-2)m-j {m - 1) = t',n-2{{'r - 1) {m - 1) - (r - 2) TO + j) = C-2(™ + 
j — r + 1) = for j = 0, . . . , r — 3. We also add a row p„ (0) = ^„^o for 
n = 0, — 2 to the table, so that the recursion ^ holds for m = 1. 
This part of the p-table, p„ (to) for < to and < n < (r — 2) (m + 1), 
is shown below for r = 4. Note that po (m) = tm-i{{r — l)m) = 1 for all 
TO > 1, and also po (0) = 1. 
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The rotated and shifted table p„ (m) and its polynomial 
extension below the staircase, for r = 4. The bold numbers 
occur also on the second subdiagonal in the table below. 



We obtained the recursion (O as a discrete integral from ^ and We 
can now take differences in recursion ([5]) and get Pn(™)~Pri-i (™) = Pn{m— 
1) - Pn-r{m - 1), or 

Pn{m) -p„(m- 1) =p„_i(m) -p„_r(™- 1) 

for all m > 1 and < n < (r — 2) (m + 1). The column pq (m) can be 
extended as a column of ones to all integers m; hence po (m) can be extended 
to the constant polynomial 1. The recursion shows by induction that the n- 
th column can be extended to a polynomial of degree n, and by interpolation 
we can assume that we have polynomials in a real variable. The extension 
of Pn (to) is again denoted by p„ (m) . The above table shows some values 
of the polynomial expansion in cursive. The expansion follows the same 
recursion, hence 

Pnix) -Pn {X - 1) =Pn-l{x) ~ Pn-r{x ~ 1) (9) 

with initial values P(r-2)m+j (to) = for j — 1, . . . , r — 2 and to > 0. These 
conditions, together with (0) = 1, determine the solution uniquely. 

Recursion ([9]) shows that (p„ (x)) is a Sheffer sequence for the same 
operator B as the sequences (s„ (a;)) in recursion ([2]). Hence pn (a;) can 
be written in terms of the same basis, the Eulerian coefficients, as s„ (x). 
However, the initial values (zeroes) for (p„ ix)) are more difficult, because 
they are not on a line with positive slope. We introduce know a Sheffer se- 
quence ((7n(a^; Q^))„>o fo'^ tl^s delta operator B that has roots on the parallel 
to the diagonal shifted by a + 1, g„(n — a — 1; a) = 0, and agrees with (p„) 
at one position left of the roots, for each n. 
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Lemma 8 For the Sheffer sequence {qn {x;a)) for B with initial values 
qo (m; a) — 1, qn{0) = Sn.o for < n < a and qn{n — a — 1; a) ~ for 
n > a, 
holds 

qn+a {n; a) = P(r^2)n-a (n) 

for all n > \a/ (r - 2)]. 



We will proof this Lemma in Subsection 13.21 
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The polynomials qn {m, 2) for r = 4 



The sequence (g„) agrees with the Euler coefficients &„ {x) 



for 



the first degrees n = 0, . 
Sheffer sequences that 



, a. It follows from the Binomial Theorem for 



qn {x\ a) 



E 

4=0 



i 



1 



x + a + 1 — nfx + a + l — i 
n — i 



X + a + 1 — i 

Corollary 9 The number of ballot paths avoiding r 1 -steps equals 

1 



(10) 



tn {m; r) 



171 — 71 — 1 

E ■ 

i=0 



1 - i 



I ^ rn + n 
i 



for m > n > 0. Furthermore, i„ (rt) = (n) = 
n > 0. 



m + 1 — i 
m — i 



n+1 
n 



for all 



Proof: <„ (m;r) = p(r„i)(„+i)_„ (n + 1) = g,„ (n + l;m - 1 - 7i). ■ 
The Corollary shows that the number of Dyck paths to (2?i, 0) avoiding 



r up steps, 



equals the number of Dyck paths to (2n, 0) 



avoiding r down steps (see formula ([4])). 
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3.1 A Conjecture for the Case r = 4. 



A Motzkin path can take horizontal unit steps in addition to the up and 
down steps of a Dyck path. Suppose a Motzkin path is "peakless", i.e., 
the pattern uu and ud does not occur in the path. Denote the number of 
peakless Motzkin paths to (n, 0) by M' (n). Starting at n = we get the 
following sequence, 1, 1, 1, 2, 4, 7, 13, 26, 52, 104, 212, 438, 910, ... for M' (n) 
(see http://www. research. att.com/~;ajas/sequences/A023431). It is easy to 

n — i\ 1 / 2i 



show that M' (n) = E[5o'^ ("2, ') m 



For the case r = 4 we conjecture that p„(0) = ( — 1)" M' {n — 3) for all 
71 > 3. That would imply E„>oP" (0) ^" = + t - ^ {1 + tf + 4t^^ /2, 
and therefore 
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n>0 ^ 

For example, the coefficient of in this power series equals 

(x - 3) {x^ + 24:x^ + 247x4 + 426x3 - 38x2 _ 2340^; + 6720) /7!, 

which in turn equals pr (x), as can be checked using the table for pr (m) 



3.2 Proof of Lemma [8] 

Because of recursion ^ we obtain the operator identity 

which holds for (p„ (x)) and (q„ (x;a)), and shows that both polynomials 
enumerate lattice paths with steps (0, 1} , (1, 1) , (2, 1) , . . . , (r — 1, 1) (above 
the respective boundaries). The number of such paths reaching (n, m) can 
also be seen as compositions of m into n terms taken from {0, 1, . . . , r — 1}. 
In the case of p„ (m) the terms ai, . . . , a„ also have to respect the boundary, 

k 

which means that Ui < (r — 2) fc for all fc = 1, . . . n. For (/„ (m; a) we get 

i=l 

k 

for the same reason that ^ bi < k + a — 1. Such restricted compositions 

1=1 

have the following nice property. 

Lemma 10 Let c G Ni, a G Nq, and n be a natural number such that 
n — a > 0. Let be the number of compositions of cn — a into n parts 

k 

from [0, c + 1] such that ai + 02 + . . . + a„ — cn — a, and ai < ck for 

i=l 

all k — 1, ... ,n — I. Let be the number of compositions of n + a into n 
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k 

parts from [0, c + 1] such that bi + b2 + . ■ ■ + bn — n + a, and hi < k + a 
for allk=l,...,n~l. Then ^ Q"^. 

k 

Proof: Suppose, 61 + ^2 + • • • + = n + a, and ^ 6^ < A: + a. Define 
ai = c+ 1 — fcn+i-i for i — 1, . . . ,n. Note that G [0, c+ 1], and n~k + a > 

n — k k 

J2bi=^n + a~J2 bn+i-i- Hence 

i=l 1=1 

fe fc 

= (c + — fori+l-i < (c + l)fc — rt + (?1 — fc) = c/c 

and 

n n 

^ fli = (c + 1)71 - ^ fcri+i-j = cn - a. 



We apply this Lemma with c = r — 2 to obtain p(j,_2)n-Q (") ~ Qn+a {'n). 
3.3 Abelization 

Let (6„) be the basic sequence for some arbitrary delta operator B, i.e., 
Bbn = bfi-i and 6„ (0) = (5o,„. Every basic sequence is also a sequence of 
binomial type, which means that X]n>o ^" (^) ~ where (3{t) = 

t + a2t^ + . . . is a formal power series. The compositional inverse of /3 {t) 
is the power series that represents B, 

B = {D)^D + b2D^ + ... 

where D — d/dx is the x-derivative. The Abelization of (&„) (by a G K) 
is the basic sequence ( f bn (x + an)] for the delta operator E^'^B 

yx+an / n>0 

(see [4]). Note that with any Sheffer sequence (s„) for B the sequence 
{sn (x + c — an)) is a Sheffer sequence for E'^B. Hence 

(( ''^x+"°" ^"('^ + '^))) >o ^ Sheffer sequence for E'^E'^'B = B again. 
Choosing c = a = 1 shows ([3]) . 

Sheffer sequences and the basic sequence for the same delta operator 
are connected by the Binomial Theorem for Sheffer sequences, 

n 

s„{y + x) Si (y) bn-i (x) . 

i=0 
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Applying this Theorem to ^ x+an ^" ~^ an)^ and (6„ (x + c + an)) shows 
that 

n 

(2/ + a; + c + an) = 6i (y + c + ai) ^ ^^n-i (x + a (n - i)) . 

^ — ' X + a in — i) 

i—O 

Choosing xasx + a + l — an, a = 1, c = 0, and y — —a — 1 gives 
K (x) = X;Lo biii-a-l) (a; + a + 1 - i).This is not quite 

what we have in pUj) : there the summation stops at a. This effect in ((TU)) 
is due to the 'initial values' 6^ (i — a — 1) which are for i > a. 

The generating function of a Sheffer sequence (s„) for i? is of the form 
(t) e-'^^*), where (t) = E„>o (0) ^ (x) ^^±;^6„ (x + c) = 
6„ (a; + c) - ;^6n (a; + c) then 

y + C) r - ^|.g(.+c)/3(t) ^ ^{x+c)m 

^-^ X + c X + c at 

n>l 

and 

J2 (a:) r = e(^+=)'3(t) _ (^^^ 

n>0 

If c = a = 1 and e^(*) - (l + t + ■ • • + T"!) , then /?' (t) = '"(^iT^f+ff^f , 
and therefore 



(1 -r)' 
(1 - 1)"'- 

the generating function of the number of ballot paths to (n, m), avoiding 



E ^» ("^) = J+2 i^t"^ (1 - + (1 - ^'^) (1 - 2i)) , (11) 

n>0 ^ ' 



4 Euler Coefficients 

The coefficients of the polynomial 



fl+t + i2 + ...+r-l) 



were considered by Euler in [2J, where he gives the following recurrence: 
n\ -J-^ f n \ fk ~ i 
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To calculate the Euler coefficients in terms of only binomial coefficients, 
we rewrite the polynomial as follows: 



{i + t + f + --- + r-y 



Thus we have proven 



1-f 



i=0 



1 - t , 

j 



Lfe/rJ 



3>0 



n\ f n + k — ri — 1 
k — ri 



(12) 



Note that this identity implies lim^^c 



fc- 1 



Combi- 



natorially, these are all {!,—>} paths avoiding — This problem occurs 
in Will's generating functionology 6 , Section 4.12. Note that identity p2p 

/n + — 1^ 

fcj„ " V k 



implies limr^c 



X 


1 


8 


36 


120 


322 


728 


1428 


2472 


3823 


7 


1 


7 


28 


84 


203 


413 


728 


1128 


1554 


6 


1 


6 


21 


56 


120 


216 


336 


456 


546 


5 


1 


5 


15 


35 


65 


101 


135 


155 


155 


4 


1 


4 


10 


20 


31 


40 


31 


20 


10 


3 


1 


3 


6 


10 


12 


12 


10 


6 


3 


2 


1 


2 


3 


4 


3 


2 


1 








1 


1 


1 


1 


1 




















1 


























n : 





1 


2 


3 


4 


5 


6 


7 


8 



A table of Euler coefficients 



for r 



We now show some properties about Euler Coefficients similar to the 
basic properties of binomial coefficients. 

1. For binomial coefficients, this property is usually called Pascal's Iden- 
tity: 

\ '■-1 / 1 



E 

i=0 



k — i 
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Proof 

r-l 

= ^f(l + i + --- + t'-i)"-i 



2. The table of Euler CoefRcients is symmetric similar to Pascal's Tri- 
angle: 

' n\ I n 

Proof We proceed by induction on n, fixing r > 2. For n = 2 we have 
the well known symmetry for binomial coefficients. Suppose true for 
some I > 2. From the above recurrence we have 

§ ((^ + -l)-k-i)^^ {{1 + l)(r - 1) - fc 
and the induction follows. ■ 

3. This property is similar to Vandermondt Convolution for binomial 
coefficients: 

+ m\ \ - / n\ / m 

k ) ^ \k~ij 

It follows because the Euler coefficients are of binomial type [J] . 

4. Here we have an identity that is trivial for binomial coefficients, i.e. 
r — 2, and gives and identity for the Catalan numbers as r — > oo. 

1 /n+i\ / n\ sr-^ . / n 



E 



n+lV n J \n I ^-^ \n — i — 1 

\ / r \ / r ^—i \ 

— -[- 1 / X ^ 1 \ 

Proof Let s„(x) = ) , as in (l3]). The binomial 

x + 1 \ n J ^ 

n 

theorem for Sheffer sequences states that s„(x + j/) — Si{y)^n-i{x). 

i=0 
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Let X = n, y = and noting that s„(0) = (1 — n) | ) = 1 — n for 

n 

< n < r and otherwise, we have s„(n) = ^ Si(0)6„_j(n), hence 

i=0 



r-2 



„ r-^ \n — I — 1 

1=1 

We have already noted that [^] _^ fn + k ^ j a,s r — > oo, so 
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